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We study the conserved charge fluctuations, as quantified by the corresponding susceptibilities, 
in strongly interacting matter as motived by the quark-gluon plasma. Using the gauge-gravity 
correspondence approach, we study the patterns of conserved charge fluctuations in two types of 
holographic models for QCD, the D4/D8 and the D3/D7 models. We compute and compare the 
quark number susceptibilities in both models and flnd an interesting common feature of the two: at 
very strong coupling higher order susceptibilities are suppressed and the conserved charge fluctua- 
tions become purely Guassian. In light of the state-of-the-art lattice QCD results we also discuss 
what we can learn from these susceptibilities about the underlying degrees of freedom in the 1 ~ 2Tc 
f^ quark-gluon plasma and examine the viability of different ideas such as holography, quasi-particles, 

fv^ as well as bound states. From analyzes of second order susceptibilities we conclude that the bound 

, states exist and are important in the 1 ~ 2Tc region. We further construct and make predictions 

C~) , for several ratios of fourth-order susceptibilities that can sensitively reveal such bound states. 
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I. INTRODUCTION 



< 

•T^ The many body systems in the strong interaction sector of the Standard Model, described by the Quantum Chro- 

O I modynamics (QCD), show very rich phase structures and provide fascinating examples of strongly interacting matter 

1^ [HE]- Amongst others a color-deconfined and chiral-symmetric phase at high enough temperature (and low bary- 

7^, onic density), known as the quark-gluon plasma (QGP), is predicted from QCD and has been created and probed 

in rclativistic heavy ion collision experiments at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron 

T— I Collider (LHC). Both the first-principle lattice QCD simulations of the QGP and the experimental measurements of 

^ its properties have consistently shown that the QGP in the temperature region T ^ a few Tc (the transition tem- 

vN perature ~ 160MeV in QCD) is still a very strongly interacting system rather than an asymptotically free plasma 

[31 3] and may be dominated by emergent degrees of freedom such as the chromo-magnetic monopoles [MZ] ■ While 

the transition to QGP at high T is a crossover, it is expected that there will be a first order phase transition at high 

baryonic density and low temperature together with a critical end point (CEP) at certain critical temperature and 

']^ density that separates the crossover region from the first order transition. There are intensive efforts both from lattice 

^rr QCD and from Beam Energy Scan (BES) experiments at RHIC looking for the critical end point on the QCD phase 

I diagram [8- 14J . To study the strongly interacting QCD matter in these relevant temperature and density regions is 

H^' theoretically challenging but of great interest and importance. 

A very useful class of observables for studying such strongly interacting matter includes quantities describing the 
fluctuations and correlations of conserved charges of the underlying microscopic theories. In the context of QCD 
;_( (counting the three light flavors), the conserved charges include the baryon number B, isospin /, strangeness S, as 

Cd well as the electric charge Q. All these charges are carried by the quarks (the fundamental fermions) or their various 
combinations (e.g. mesons/baryons as bound states of quarks), and the patterns of their fluctuations/correlations 
provide sensitive probes to the actual degrees of freedom that carry such charges |15U16| . Such patterns are particularly 
important imprints of potentially emergent degrees of freedom when the system is strongly interacting: a famous 
example from condensed matter physics is the fractional quantum hall state with emergent quasi-particles carrying 
fractional electron charges that were experimentally first identified through current fluctuation measurements. For 
simplicity we focus in this paper mostly on the baryon number fluctuations which is also simply related to the 
quark number fluctuations as each quark/anti-quark carries ±1/Nc baryonic number and the corresponding baryonic 
chemical potential and quark chemical potential are related as fj,B = NcfJ-q. Let us consider the thermodynamics 
of QCD-likc matter with Nf flavors of fundamental quarks carrying the baryonic charge, described by the pressure 
P(T, /i) as a function of temperature T and quark chemical potential /i. These fluctuations can be quantified by the 
susceptibilities, defined through the Taylor expansion coefficients of pressure over the chemical potential [T7H2T| : 






The above dimensionless susceptibilities are closely related to charge densities fluctuations and correlations: Xf ^< 
B > (measuring density itself), X| '^< B'^ > (measuring quadratic fluctuations), etc. When there are multiple 
charges, the susceptibilities can be easily generalized to all charges, and there will also be cross-charge (ofl'-diagonal) 

susceptibilities that describe the density correlations, e.g. Xbs — o( , /T)d( /t) '^< B S >. Of particular interest 
is the zero density limit of these susceptibilities which can be directly evaluated by lattice QCD simulations, and we 
also introduce the following zero density susceptibilities (per quark flavor for convenient comparison with lattice) 

,m.'^^^^ (2) 

Note in the above by symmetry all the odd susceptibilities vanish at /i — > and only n = 2,4,6, .. terms survive. 
These susceptibilities are also closely connected with many observables in the heavy ion collision experiments, such as 
the charge fluctuations and correlations measurements, the search for critical end point, and the freeze-out conditions, 
etc [HI [^2^EB] . Let us just mention one such example: the cumulants of net baryonic charge fluctuations proposed 
for the search of CEP. Let a be the variance and k the kurtosis of the event-by-event net baryon number distribution, 
then for a thermal system one has the following relation 

^^ ^ ^Kt>) ^'^ 

It is therefore also interesting to study such ratios of susceptibilities. 

One significant challenge to study these susceptibilities for the quark-gluon plasma in the physically interesting 
regime, as already mentioned, is that the system is strongly interacting. Indeed the state-of-the-art lattice QCD 
results [ini m] show rather distinctive patterns in the 1 ~ 2Tc region that can not be understood via two often-used 
simple benchmarks, namely the free hadronic resonance gases at low temperature and the Stefan-Boltzmann gas of 
quarks and gluons at asymptotically high temperature (see e.g. discussions in [27j). To this end, the holography 
provides a very useful approach for understanding the behaviors of the conserved charge fluctuations in QCD-like 
strongly interacting matter. As is well known, the holographic models in the past have offered many insights for 
understanding the properties of the quark-gluon plasma (see e.g. the comprehensive review |28j and references 
therein). There were though only very few studies of these susceptibilities in holographic models using one or another 
setup [371 [Ml ISO] ■ In the main part of this paper, we will use the gauge-gravity correspondence approach to study 
the patterns of conserved charge fluctuations in two major types of holographic models for QCD, the D4/D8 and the 
D3/D7 models in the Section II and III respectively. We will compute and compare the quark number susceptibilities 
at various temperature and density in both models and also study the related quantities such as the ratios of them. 
We will analyze and compare the asymptotic behavior at large coupling and large temperature in both models and 
find interesting common features and differences. Another part of the motivation for this work, is to closely examine 
the implications of the state-of-the-art lattice QCD results for the susceptibilities, which will be done in the Section 
IV. We will evaluate the viability of different ideas such as holography, quasi-particles, as well as bound states, and 
discuss what we can learn from them about the underlying degrees of freedom in the 1 ~ 2Tc quark-gluon plasma. 

II. SUSCEPTIBILITIES FROM HOLOGRAPHY: D4/D8 MODEL 

In this section we will evaluate the susceptibilities from one holographic model for QCD based on the D4/D8 branes 
configuration with one compactified dimension, known as the Sakai-Sugimoto model |31j . In this model, the flavor 
dynamics is given by Nf D8 — D8 flavor branes in the background flclds (the gravity, dilation, and Ramond-Ramond 
fields) generated by Nc DA branes in the "probe" limit Nj « N^. The Sakai-Sugimoto model at zero temperature 
has been shown to reproduce many results in qualitative or semi-quantitative agreement with QCD ^31] . The model 
has also been extended to finite temperature [32] and finite baryonic density [33] , with its phase diagram qualitatively 
mimicking that in QCD. We refer the readers to the above references for detailed descriptions of the model. 

Let us now discuss the different phases at finite T and fi in the Sakai-Sugimoto model. At flnite temperature 
there are two possible background geometries (the zero temperature one and the black hole one) and two important 
scales involved, the Kaluza-Klein mass scale Mkk (from the compactified dimension) and the temperature scale T 
(ultimately associated with the black hole horizon). At low temperature the zero temperature geometry has the 
smallest action and is thermodynamically preferred: this is the confined phase. When temperature is high enough 
T > Tc = Mkk/{2tt) the black hole geometry will "win" and the system is in the deconfined phase. By turning on 
a gauge field in the bulk one induces global conserved (baryonic) charge density (and the related chemical potential) 
in the boundary field theory and mimics the finite density. This does not change the confined/deconfined phase 
distinction determined by temperature alone. In the confined phase, however, when the chemical potential is high 



enough fj, > fie = XMKKf^^n (with A = g'^Nc the 't Hooft couphng), there is a new phase (the cold dense phase) 
that has nonzero charge density and resembles the proposed "quarkyonic phase" [34] . In the "vacuum phase" T < Tc 
and n < He thermodynamic quantities are independent of T, ^ and therefore all susceptibilities vanish. We will focus 
on the "quark-gluon plasma (QGP) phase" at T > Tc and the "cold dense phase" at T < Tc and fi > ^c and obtain 
the susceptibilities in each phase respectively. The phase diagram is demonstrated in Figjl](left panel). 
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FIG. 1: The phase diagrams of the D4/D8 model (left) and the D3/D7 model (right, with the inset plot at the top right corner 
showing the details close to Tc). 



For convenience and simplicity of the presentation, we hereby define some dimensionless variables as following: 
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In the above, the P is pressure, and d is the quark number density at a given chemical potential /i ( — noting that the 
baryon number density would then be pB = d/Nc and the baryonic chemical potential fiB = Nch). In addition there 
are two parameters R = {ngsNcl'^y^^ with gs, Is being string coupling and length, and Ut = (47r/3)^i?^r^ the horizon 
radial coordinate in the black hole geometry. More detailed definitions and computations in the Sakai-Sugimoto model 
specifically related to our discussions can be found in e.g. [27,. It should be pointed out that for the D4/D8 model 
we will closely follow the approach of [27] in which results for susceptibilities at /i = in QGP phase and bX fi — jic 
in cold dense phase were reported. Here we will compute the susceptibilities generally for any given T and /i. 

A. Susceptibilities in the QGP phase 

For the QGP phase of the D4/D8 model [32l [33], the pressure has been computed at given temperature T and 
density d to be (with 2^1 the hypergeometric function): 
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where the density d(T', fi) can be determined from the foUowing relation: 
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With the chain rule, §£ — „_(„- and generally ^~„J^ — — — S ,t.~ — , we can compute the susceptibilities to 

<^M a^/ad *^M a^/od 

arbitrary order defined as: 

fi^(p/T*) ~ 8"P 

X«(T, fi) = J ' ' = (33"-«24-"^-i) iV^TV/A^-" r"-4-^ (9) 



A few explicit examples, for the first to fourth order, are given here: 
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Let us first discuss the results at /i = for which only even-order susceptibilities survive. In this case we have 
simple analytic results: 



XI{T) = -^^ '- = ?„A3-" (^^ j (14) 

with ^„ being numerical constants, e.g. ^2 ~ 0.071, ^4 « 2.2, and ^g ~ —158 for Nc = Nf = 3. These results were 
previously obtained in |27j . 

We then examine the patterns of X^^ at nonzero fi. In Fig[2| wc plot ^^^1 23456 ^^ ^ function of T > Tc for a 
variety of values /i. The asymptotic behavior of the first andsecond order susceptibilities could be understood as 
follows: analytically one can show that at high T limit one has the leading contribution P ~ fjl^T^ and therefore 
Xf ^ ub/T^ (x fiT'^ while X| ex /x" T. All the higher order susceptibilities also show interesting patterns and in 
particular all are suppressed with increasing temperature. We emphasize that these results are useful and provide 
direct information on the fluctuation patterns in a strongly interacting matter at regimes with /i comparable or lager 
than T where the Taylor expansion via small fi/T is not useful. 

Let us now discuss the ratio ^ (in unit of T^ /fj/^ = (27/(2A))^) which is a quantity of particular interest for the 
CEP search for a wide range of T and /x values in the QGP phase: see FigJ3| As one can see from the plots, the 
general trend is that this ratio decreases when getting to larger T and /i values and is always positive. If one takes the 
phase boundary in this model (i.e. the T — T^ line) as an analog to the freeze-out boundary in heavy ion collisions, 

then walking from small ^ toward large /x one sees a reducing value for —^ which shows a qualitatively similar trend 

^2 
to model computations of QCD and to preliminary STAR measurements. 

B. Susceptibilities in the cold dense phase 

In this subsection we study the susceptibilities in the cold dense phase of D4/D8 model. For this phase, the pressure 
and chemical potential are given by the following equations: 
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FIG. 2: -^^^=1, 2, 3, 4, 5, 6 in unit of {Tc/^c)" ^ versus T /Tc at various values of /i in the QGP phase of the D4/D8 model. The 
inset at top left Conner of X| panel shows the details close to Tc . 
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FIG. 3: -r4- versus T at various values of ^ (left) and versus p at various values of T (right) in the QGP phase of the D4/D8 
model. 
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and both quantities are independent of the temperature. We then compute the susceptibiUties with the same definition 

as in En. (9). To do that, we introduce the following notations: u*^*^ = i# and P'-'^ — ^M^. By chain rules we have 

,p(i) 



that d = 3^ ~ ^WT' ^^"-^ similarly we can get the following relations: 



3F(i) = Jl^'^d, 3P(2) = M^'^rf + M^'^ , 3P(3) ^ ^(3)d + 2/7(2) ^ 3p(4) ^ ^(4)^+3~(3) 
With these relations we can then easily derive 
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where the p"^' are given by the integrals: 
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FIG. 4: 



X^^j 2 3 4 5 6 i'^ unit of (Tc/fic)" ^ versus ^ at T = Tc in the cold dense phase of the D4/D8 model. 



In Figll^we show these susceptibilities X'^^^ 23456 ^^ ^ function of fi. While the first three susceptibilities appear 
to monofonically grow with chemical potential, the higher order ones get suppressed at large /i. The asymptotic 
behavior could be qualitatively understood analytically through the following: at large density the density and 
chemical potential is roughly related as d ~ p^^ and thus X' ^ p^^'^'" so we see the Xf 23 to be positive and 



growing, X| ^ ^ ^/^ and Xq ^ ^ ^/^ to be positive but decreasing, while X| ~ 
decreasing in magnitude. 

For the ratio ^ (in unit of T|?/^c = (27/(2A))^), we have the following formulae: 
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The results are shown as a function of fi (noting that there is no T-dependence in this cold dense pahse) in Figls] 
For comparison we also show the same ratio in the QGP phase as a function oi ij, a,t T — Tc- While both curves show 
similar magnitude and similar decreasing trend toward large /z, there is a quantitative discontinuity across the phase 
boundary. For example as labeled in the figure, this ratio at {T^, i^c) is 0.133 in QGP phase and 0.123 in cold dense 
phase though the difference is rather small. 
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FIG. 5: ^ versus chemical potential fi in the cold dense phase of the D4/D8 model and comparison with the same quantity 
aXT = Tc from the QGP phase of the D4/D8 model. 



III. SUSCEPTIBILITIES FROM HOLOGRAPHY: D3/D7 MODEL 



In this section we will evaluate the susceptibilities from another holographic model for QCD based on the D3/D7 
branes configuration [321 1351 |Mj • In this model, the background geometry is provided by the Nc D3 branes representing 
the gauge field dynamics while the fundamental matter is represented by Nf << Nc ("probe limit") flavors of D7 
branes embedded in the background geometry. These flavor branes are placed at a distance of Mq in the holographic 
dimension from the D3 branes and the flavor chemical potential /i is implemented via a U{1) gauge field in the bulk. 
Different embeddings of the flavor branes correspond to different phases of the model. On the finite T and /i phase 
diagram as shown in FiglT] (right panel), there are two distinctive phases: one at low T and /i corresponding to a 
"Minkowski embedding" with zero density, while the other at high T and /i corresponding to a "black hole embedding" 
with nonzero density, with the two separated by a phase transition line. In this study we focus on the susceptibilities 
in the nonzero density phase that has nontrivial dependence on the chemical potential fi. We will closely follow the 
setup and conventions in [35] where one can also find a detailed treatment of the D3/D7 model. For later convenience 
we also introduce the following definitions and conventions from this model: 
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with parameters ls,gs being the string length and couphng, A = g'^Nc = 2TTgsNc the 't Hooft couphng, uq the horizon 
parameter and T^v the tension of the D7 branes (see more details in [35]). Here we study the finite density phase 
above the transition line corresponding to the black hole embedding. 
The susceptibilities we would like to calculate are defined as 
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with the relation between the chemical potential /i and density d determined by 
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In the above the factors f — 1 — p^^ and f — 1 + p^'^ are from metric (with p the radial coordinate). The function 
x{p) describes the embedding of the D7 branes in the geometric background (with x — ^x/ dp) and is solved from the 
equation of motion belowf35j: 
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One can determine the phase boundary {T,p)c by solving the Eqs.( 25|26 l in the limit of vanishing density d — > as 
shown in Figfl] (right panel). Here we define Tc and pc as the crosspoint of the phase boundary line and T,p axes, 
respectively. One can get pc = Mq and T^, = 0.764 x —j=- as in |55] . 

Our strategy to compute X'^ in Eq.(24| is to start with density d as d'^P/dp'^ ^ d^^^d/dp"^^. To do that let us 
first introduce p^"^^ = ^-^ and the integrals li defined as 
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With these integrals, the Eq.(25) can be simply written as p — 2\f2 x (0.764T) d/i. One can further obtain the 
following results by differentiations: 
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Finally we can get the following formulae for the susceptibilities : 
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These susceptibilities can be analytically obtained in the very high temperature limit T » Mq and vanishing 
density d — > 0, for which one simply has x ^ 0- In tbis limit we have 
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and can get the following results for susceptibilities: 
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FIG. 6: -'f'^i 2 3 4 5 6 i'^ urii^ of {Tc/ fic)^ ^ versus T at /^ = (thin solid blue), ^ = 0.3/^c (dashed red), fj, — fic (thick solid 
magenta), and /i — 2^c (dotted black) from the D3/D7 model, respectively. 




FIG. 7: -^ versus T at various values of n (left, from bottom up with fJ,—5, 4, 3, 2.5, 2, 1.5, 1, 0.7, 0.3, 0.1 ) and versus fi at 
various values of T (right, from bottom up with r=0.3, 0.5, 0.9, 1, 1.5, 2, 2.5, 3, 4, 5) in the D3/D7 model. 



10 

Generally for any given temperature and density one can solve the equation of motion for x(p) numerically and 
then compute the susceptibilities as in the above formulae. In FigJ6^ we plot X'^^-^ 23456 ^^ ^ function of T for a 
variety of values /x. A first distinction from the D4/D8 model is that the susceptibilities diverge when approaching 
the phase boundary. A second interesting pattern is that all the susceptibilities appear to approach certain constant 
values toward large temperature in contrast to the results from D4/D8 model that show strong T-dependence as 
T — > 00. We will discuss this distinction more in the next Section. Finally we notice that the odd susceptibilities for 
/i = vanish as required by symmetry. 

Finally we discuss the ratio ^ (in unit of T^/Hc = 0.764^ • 4/A) which is shown for a wide range of T and /n values 
in the D3/D7 phase: see FigJT] A most distinctive feature here, as compared with the D4/D8 results, is that the ratio 
diverges when approaching the phase transition line from either T or n directions. Also different from the D4/D8 
case is that (for non-divergent regime T > Tc and/or fi > fj,c) the ratio here increases with temperature, though it 
decreases with /i in both models. 

IV. DISCUSSIONS IN LIGHT OF CURRENT LATTICE QCD RESULTS 

Lattice QCD simulations provide the ultimate answers for the precise values of these susceptibilities in real QCD 
matter (albeit limited to the very small density region due to the sign problem). To understand the implications of 
these results, however, could be tricky particularly in the strongly interacting regime. Current lattice QCD results 
(mostly for the second order susceptibilities) from both Wuppertal-Budapest group [2Qi and BNL-Bielefeld group [21] 
are at almost physical masses, with high precision, and in agreement with each other. In light of these lattice input, 
we discuss in this section what we can learn from various models about susceptibilities and whether these models 
provide viable and consistent descriptions for one or more aspects of the lattice QCD data. 

A. What can we learn from holographic models? 

With the susceptibilities being computed from both D4/D8 and D3/D7 models, it is tempting to compare the two 
and find common as well as different features. Let us focus on the susceptibilities Xn=2.4 6....(^) ^^ ^^™ density. 

One very interesting common feature of the two models is that all higher order susceptibilities n > 2 are suppressed 
by inverse powers of strong coupling A. Namely, we have 

X^cxA^-" {Di/D8) (35) 

X^cxA^-t {D3/D7) (36) 

Clearly in the infinity coupling limit A — >■ cx) only the second order xl survives in both models, which seems to suggest 
that the fluctuations of conserved charges become Guassian at very strong coupling. This feature might be universal, 
and if so would be the most interesting lesson we learn from holographic models. 

We can compare the temperature dependence of Xn '^^ the two models, as shown in FigjSlfor n = 2,4,6. We also 
show the lattice results for X2 (light quark number susceptibility) [201 El] . These plots are obtained with coupling 
parameter A = 9 for both models and one may note that such a choice of A may not be the optimal one for other 
physics considerations in these models and also the optimal choice of A may not be necessarily the same for the two 
models. In fact the higher order susceptibilities are rather sensitive to the choice of A. While it is hard to consider 
these results as quantitative explanation of the actual QGP susceptibilities, the qualitative trends in such holographic 
models appear rather interesting and informative. 

Finally we would like to compare the high temperature limit of the Xn for both models: 

/ rp \ 3 — n 

xUT^^)^i^-J {DA/D8) (37) 



X«(T^oo)«(^-j {D3/D7) (38) 

There are interesting differences between the two models: the D3/D7 model results show constant asymptotic values 
at high temperature for all orders (as also evident from the plots in Figjs]), which may be understandable as there 
is essentially only one scale for the model in the limit T » Mq] the D4/D8 model results however show strong 
temperature dependence which may be due to the two scales in the model even at high T limit (with more and more 
high KK models involved), and furthermore at high T all higher order n > 2 susceptibilities are suppressed and the 
fluctuation becomes Guassian. 
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FIG. 8: Comparison of Xn=2 4 eC^) from the D4/D8 model (red curves) and from the D3/D7 model (blue curves). Also shown 
for X2 ''■re lattice QCD results (black diamonds). 



B. Can quasi-particle models describe the conserved charge fluctuations? 



There have also been a lot of studies in the Hterature trying to understand these susceptibilities in terms of quasi- 
particle models [HI [37U42J . Such models are based on the idea that the (complicated) effects from interactions 
in the many body setting may be approximately accounted for by a simple medium modification of single particle 
(the quasi-particle) properties, e.g. a temperature and density dependent mass. In the context of QGP, the quasi- 
particle models assume the plasma to be a gas of independent quasi-quarks and quasi-gluons whose masses change 
with temperature and density. At very high temperature, this may indeed be a very good approximation and the 
quasi-particle properties are computable from perturbation theories, as demonstrated e.g. in the hard thermal loop 
approach |43fl45| . In the regime T ~ Kqcd, however, there is no a priori reason for a simple quasi-particle picture 
to be a good approximation. The susceptibilities are sensitive to the properties of the conserved charge carriers (e.g. 
quasi-quarks), and one may expect them to provide a useful test of the quasi-particle models. Indeed lattice results 
for the susceptibilities in the 1 — 2Tc region show rather distinctive patterns that can not be easily explained by e.g 
perturbative calculations. In an unconstrained quasi-particle model, though, one may still "adjust" the dependence 
of the quasi-quark mass Mq{T,fi) — noting that d^AI/dfi^ could provide additional contributions to susceptibilities. 
In fact within such a quasi-particle model one can generally establish a one-to-one correspondence between xt ^'^d 
A/(T, 0), xl and dAI/dfi{T, 0), xt and d'^AI/dfi'^{T, 0), etc, as pointed out in ^15j. Recently, high precision lattice data 
for susceptibilities of three conserved charges B, Q, S (baryon number, electric charge, strangeness) in QGP became 
available [201 El] , and with these three quantities a much more stringent consistency test is in order for quasi-particle 
models which have only two adjustable quasi-quark masses for the light quarks and for the strange quark. 

Let us consider such a quasi-particle model with in-medium masses M; (T, /i) for u,d quarks and Ms (T, fi) for s quark. 
To compute the fluctuations oi B,Q, S we start with their densities by summing over all quasi-quark contributions: 



d'^p 



d'^p 



1 



1 



ns = 2.3-(-l) 



(27r)3 



1 1 

1 1 



A^/^^+Mi-^^s)/T j^ ^ g(^p2+M2+^,)/T ^ ^ 



(39) 
(40) 
(41) 



with qf is the quark electric charge with g„ = 2/3 and qd,s — —1/3. Note that the flavor- wise chemical potentials are 
related to the conserved charge chemical potentials via: /j.^ = /is/3-(-2/iQ/3, ^a = ^J'B|i—^^Q|'^^ Ms = t^s/'^—IJ-Qli—l^s- 
The susceptibilities can then be computed via: 



X^{T) 



driB 



dfii 



all /i's-i-O '^'"Q 



all p's^O ^^S 



(42) 



all fi' s^O 
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Note that by charge conjugation symmetry the quasi-particle masses Mf{T,Hf) is an even function of fif i.e. 
dM/d^\^^o -^ 0. Therefore one sees in the present approach that only Mi{T) (u,d quasi-quark masses) and Ms{T) 
(s quasi-quark mass) will appear in xf ^^^ xi while the Ms{T) alone will appear in xf . This situation provides an 
opportunity to check the consistency of the quasi-particle model: one can extract the M/(T) and Ms{T) from lattice 
data for X2 ^^^ X2 ^ ^^'^ alternatively one can also extract the MgiT) from lattice data for xf , and then confront the 
two ways of extracting the same quantity Mg (T) . Using the lattice data from [20] we have done this exercise and the 
results are compared in Fig[9j While the two agree for temperatures greater than ^ 2.5Tc, one clear sees a divergence 
of the two toward lower T and the discrepancy becomes huge when T ^ Tc- 



5^ 




150 200 250 300 350 400 
T (MeV) 

FIG. 9: Comparison of strange quasi-quark mass Ms{T)/T extracted from X2 & Xs (r^d solid line with box symbols) with the 
same mass extracted from xf (blue dotted line with diamond symbols). 

From this analysis using the state-of-the-art lattice QCD data, we conclude that simple quasi-particle models can 
not provide a consistent description of the conserved charge fluctuations in the 1 — 2Tc regime. There are of course 
varied ways of implementing quasi-particle models (e.g. starting from pressure rather than from conserved charge 
densities and including T, ^ dependent Bag terms, adding Polyakov loops, etc) which may yield somewhat different 
results. But the main message is unlikely to change, namely there are very strong cross-flavor correlations (described 
by off-diagonal susceptibilities) in this regime that are generally missed by quasi-particle models. 



C. Do off-diagonal susceptibilities imply bound states above Tel 



Finally let us discuss the off-diagonal susceptibilities, e.g. Xnd, Xus which provide direct information on the density- 
density correlations between particle species carrying different flavors. In quasi-particle models such cross-flavor 
correlations are missing, while in holographic models (and generically in large Nc limit) such off-diagonal susceptibili- 
ties are 1/Nc suppressed as compared with the diagonal ones as recently shown by Casalderrey-Solana and Mateos in 
[30j . They are however very sensitive to the underlying degrees of freedom, in particular possible composite objects 
(such as mesonic states) that naturally carry multiple flavors. Since the QGP in the 1 — 2Tc regime is a strongly 
interacting system and the interactions between quarks/anti-quarks are still strong, it is plausible to have bound 
states (both meson and baryon like ones, and possibly colored bound states) above T^ [46l [47]. The off-diagonal 
susceptibilities thus can be a good test of the existence of such bound states. In fact observables along this line were 
proposed quite some time ago: one example is the B-S correlation Cbs = — 3 < BS > / < 5^ > [T0| which is a sensi- 
tive probe of the strangeness carrying states; another example is the high-order B-I (I being the isospin) correlations 
< B'^P > / < B^ >< P > and < B^P > / < B^ >< P > [15] which are particularly sensitive to existence of high 
baryon charge objects such as the survived baryonic bound states above T^. Recent high precision lattice data e.g. 
for the Cbs indeed show significant departure from simple independent (quasi-)quark picture and the nice analysis 
in [48 appears to suggest that bound states are indispensable for explaining such off-diagonal susceptibilities. 
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FIG. 10: (left) lattice QCD results for Xv.a\ (middle) extracted temperature dependent mass M^{T) for kaon states; (right) 
contributions of these kaon states to higher order off-diagonal susceptibilities x-u-uss- 



Mesonic states 



Let us first focus on the lattice results for Xus [10] as shown in Fig 10 (left). Note that for a gas of quarks at high 
T one would expect Xus — ^ while in the hadronic phase T < Tc ~ 150MeV the Xus can be explained by a hadronic 
resonance gas model. There is strong negative contribution even above Tc ~ IbOMeV and continuing toward ~ 2Tc 
— note also that if one simply continues the vacuum hadrons' contribution to T > Tc the curve would then continue 
going more and more negative rather than bending back. A possible explanation of this behavior could be the gradual 
melting of hadronic bound states with increasing temperature [13 \M1 HI] • Among the hadronic states contributing 
to Xus, the two charged kaon states (K^) are the lightest and most robust, and indeed they contribute negatively 
to Xus- For simplicity we assume these states dominate the bound states' contribution to the Xus, and examine the 
implication of the lattice data. To mimic the "melting" of these states, we assume an effective mass Af^ (T) for these 
kaon states, and their contribution to the Xus is then given by: 






27r2 



dx- 



2gV^^ + (M^/T)2 



X e 



^^x-^HM'^ITY 
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(43) 



Then from the lattice data on Xus one can extract the effective mass M'^^T^, as shown in the Fig 10 (middle). As one 
can see, the mass is essentially the physical vacuum value for T < Tj,, while starts to rapidly increase for T > T^, which 
means these kaon states are more and more strongly suppressed toward higher temperature. Of course such increasing 
mass is just an effective way of mimicking the melting of such states due to increased screening of interaction potentials 
and in reality it may not be plausible to think of the kaon-like states with ^ GeV mass. To see if such a way of 
mimicking the melting of bound states above T^ is a good approximation, we can use the higher order off-diagonal 
susceptibilities to test it. Contributions of these kaon states to the Xuuss (and Xuuss = Xusus — —Xusss = —Xuuus 



from these states) based on the present assumption and extracted M^{T) are given in the Fig 10 (right), which 
can be tested by future lattice QCD data. We nevertheless point out that there are in principle contributions from 
various other mesons (e.g. excited charged kaon states) and baryons (e.g. S, A, S states) and some of the baryons' 
contributions may become more visible in higher order susceptibilities (as S, S carry more units of u or s flavors). On 
the other hand, if these higher mass, more loosely bound states could survive above Tc, there is no reason why the 
lowest kaon states wouldn't. 

Another place where the strange mesons (mainly kaons) can play an important role is the strangeness-electric charge 
correlations. On the quark side only strange quarks/anti-quarks contribution while on the hadronic side the kaons 
will contribution, and importantly the K have their electric charges three times more than the quarks so they will be 
more and more prominent in higher and higher order S-Q correlations. For example, one can construct the following 
ratios of fourth-order susceptibilities 



C, 



22/13 
QS 






C. 



31/13 
QS 



9 



<Q53 > 



(44) 



There are normalized such that the strange quarks/anti-quarks will contribute unity. The K mesons will contribute 
3 and 9 respectively and will be more visible in the latter ratio. One can make a simple model to compute these 
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ratios. If we assume at the second order susceptibilities in QGP, the S'^ is dominated by strange quarks while the Xus 
is dominated by the K^ mesons, then we can use the extracted Ms{T) and M^{T) to make predictions for these two 
ratios: the results are shown in Fig{TT] The patterns are quite clearly showing a gradual switch from meson-dominated 
contributions to the quark-dominated contributions from low to high temperature and the switching is slower in the 
where the mesons should be more visible. At very quantitative level, the fraction from strange quarks in the 



31/13 
QS 



plots may be a bit overestimated for two reasons: first the < S^ > should not be entirely due to strange quarks; 
second there are many other higher strange mesons and baryons states which will contribute. So the actual ratios 
may be slightly higher than the simple model estimations here, but we expect the patterns to be readily there and 
the numbers to be quite quantitative. 




250 300 
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400 



FIG. 11: The ratios C|g/" 
for detailed discussions. 



(blue diamond) and Cng (red box) for the 4th order baryon-electric charge correlations: see text 



2. Baryonic states 



A more direct probe to "tell" if there is any baryonic bound state would be the direct correlations between baryon 
number and other conserved charges such as the baryon-isospin correlation proposed in [15| or the baryon-strangeness 
correlation proposed in [TS]. Note however that in the case of baryon-strangeness Cbs = — 3 < BS > / < S^ >, there 
are still strange mesons' contributions to the < 5^ > in the denominator. (Similar issue exists for the baryon-isospin 
correlation as the normalization < /^ > involves mesons' contributions.) In order to probe the purely baryonic states 
and avoid the complication with the mesons, one can use specific fourth order cumulants. For the baryon-strangeness 
correlations, one can study the following ratios: 



C 



22/13 
BS 



<BS3> 



C 



31/13 
BS 



9 



< g3gi > 
<BS^ > 



(45) 



In these ratios, only strange baryons and strange quarks contribute to both the numerator and the denominator 
and they are normalized such that the pure strange quarks will give unity for both. In Fig 12 (left and middle) we 



indicate the values of these ratios for a single component made of various baryonic states based on their B and S 
quantum numbers. It is more difficult to make a quantitative estimate for these ratios, as the observables involve a 
lot of baryonic states with varied masses and medium modifications. We can nevertheless make a simple estimation 
for the region not too much above Tc by a mixture of all baryonic states in Particle Data Book below 2GeV and 



the results are shown in the Fig 12 as the green bands. If indeed the baryonic states survive above Tc and gradually 
melts away, then these ratios should gradually switch in the 1 — 2Tc regime from the baryon-dominated values (as 
indicated by the lines for the baryonic states) to the quark-dominated values (the lowest line at unity). Furthermore 



31/13 



-,22/13 



the states with larger baryonic charges (the baryons) are enhanced more in the ratio Cgg compared in the Cgg 
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31 /13 

SO the approaching to the quark value with increasing temperature will be slower for C^g . 
be readily tested by the accurate lattice QCD data for 4th order ofF-diagonal susceptibilities. 



All these patterns can 
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FIG. 12: 



Various baryonic and quark states' contributions to the proposed 4th order baryon-strangeness Cg^ (left), 



^31/13 

(middle) and baryon-isospin correlations Cg/ (right), with the green bands indicating the value for the mixture of all barynoic 
states in Particle Data Book below 2GeV: see text for detailed discussions. 



Similarly for the baryon-isospin correlations, one can study the following ratio: 

9 < B3/1 > 



C 



31/13 
BI 



4 < BP > 



(46) 



In the above ratio, only baryons with nonzero isospin and the u,d quarks contribute to both the numerator and the 
denominator and they are normalized such that the pure u,d quarks will give unity for both. In Fig 12 (right) we 
indicate the values of the ratio for a single component made of various baryonic states based on their B and S quantum 
numbers, and again we make a simple estimation for the region not too much above T^, by a mixture of all baryonic 



states in Particle Data Book below 2GeV and the results are shown in the Fig 12 as the green band. Note for the 

22/13 

isospin it is slightly more complicated and the C^/ would be rather complex due to the opposite contributions from 
isospin partners (such as the protons and neutrons as well as the u and d quarks) . We expect this ratio will gradually 
decrease from the green band and change toward the pure quark value when temperature is increased to ^ 2Tc. 

Finally it is not difficult to generalize such ratios to the 6-th order susceptibilities that provide even more sensitive 
tests to the existence of bound states which future lattice data will be able to prove or disprove. 



SUMMARY 



In summary we have studied the conserved charge fluctuations, as quantified by the corresponding susceptibilities, 
in strongly interacting matter via a number of approaches. The holographic models motivated by QCD matter provide 
ways to overcome the challenge posed by the very strong interaction in many body setting, and we have used two 
types of such models, the D4/D8 and the D3/D7 models, to explore the patterns of conserved charge fluctuations 
in these systems. We have explicitly computed in both models the susceptibilities and certain combinations of them 
that are of physical interest, with both analytical results when possible and numerical results when needed. While 
the susceptibilities from the two models show different behaviors in several aspects, they share the common feature 
that at very strong coupling higher order susceptibilities are suppressed and the conserved charge fluctuations become 
purely Guassian. This has not been known before, and it would be very interesting to see if this feature is universally 
true also in other very strongly coupled systems. 

In light of the state-of-the-art lattice QCD results for (mostly second order) susceptibilities with distinctive patterns 
in the 1 — 2Tc region, we have also examined the viability of different models aiming to explain these data. The 
holographic models show qualitatively interesting trends that are similar to lattice results but may not provide a 
quantitative explanation. The quasi-particle models are shown to contain inconsistency in explaining all conserved 
charge fluctuations due to the strong cross-flavor correlations that are missing in such models. It has also been 
demonstrated that hadronic bound states, which survive above Tj, and gradually melt away, provide a reasonable 
description of the off-diagonal susceptibilities. More future lattice QCD data for the higher order susceptibilities and 
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particularly for the off-diagonal ones will be crucial for discriminating different models and for understanding the 
degrees of freedom in the strongly coupled quark-gluon plasma. 
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